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I present an extensible experimental design for optical continuous-variable cluster states of arbitrary size 
using four offline (vacuum) squeezers and six beamsplitters. This method has all the advantages of a temporal- 
mode encoding [Phys. Rev. Lett. 104, 250503], including finite requirements for coherence and stability even as 
the computation length increases indefinitely, with none of the difficulty of inline squeezing. The extensibility 
stems from a construction based on Gaussian projected entangled pair states (GPEPS). The potential for use of 
this design within a fully fault tolerant model is discussed. 

PACS numbers: 03.67 .Lx, 42.50.Ex 
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I. INTRODUCTION 

Quantum computation (QC) is the controlled coherent ma- 
nipulation of quantum information. Analogous to its classical 
counterpart, quantum information is usually encoded in the 
quantum state of locally addressable systems, and its manipu- 
lation is performed by inducing coherent unitary evolution of 
these systems via external laboratory equipment [ 1 1. The one- 
way model of QC J2[ replaces this need for coherent unitary 
control by a sequence of adaptive local measurements made 
on a highly entangled resource state called a cluster state Q. 
This resource acts as a universal substrate with quantum in- 
formation encoded virtually within it. With local projective 
measurements often being easier to implement than unitary 
evolution, the philosophy of the model is this: perfect the cre- 
ation of cluster states, and the rest is (comparatively) easy. 

In its continuous-variable (CV) incarnation, universal one- 
way QC [4| requires a resource state known as a continu- 
ous variable cluster state [5|, which is a multimode squeezed 
Gaussian state. In an optical setting, homodyne detection and 
photon counting — plus classical feedforward of measurement 
results — suffice to implement universal QC using CVs J6). 
Homodyne detection alone is sufficient to implement all mul- 
timode Gaussian operations 16 1, given a cluster state with a 
sufficiently connected graphPl 

There are currently four proposed methods for construct- 
ing universal CV cluster states optically [4 9-11 1. These all 
have the advantage of deterministic preparation over their op- 
tical qubit counterparts fl^HHL which rely on nondeterminis- 
tic interactions and postselection to generate the desired state. 
Each has specific advantages. 

The canonical method [4| involves single-mode squeez- 
ers and controlled-Z (Cz) gates, which are an example of a 
quantum nondemolition (QND) interaction [ 15 1. The Cz gate 
can be implemented using beamsplitters and inline squeez- 
ing (i.e., squeezing of a state other than the vacuum) |fl6l[T7l . 



1 While usage varies in the literature, I am using the convention |4 7 1 1 that 
a "cluster state" can have any associated graph. Some authors prefer to call 
this a "graph state," reserving "cluster state" for square-lattice graph states 
only. I would refer to such states as "cluster states with a square-lattice 
graph." 



which is experimentally challenging but achievable using cur- 
rent technology |18|. All Cz gates commute, and this gate 
is the natural CV generalization of its qubit counterpart (used 
to create the qubit cluster states [3|), thus making this design 
particularly amenable to theoretical analysis. 

The linear-optics method [ 1 1| eliminates the need for inline 
squeezing by replacing the Cz gates with a network of beam- 
splitters. The replacement is not one for one, however, and an 
entirely new network will generally be needed to make a clus- 
ter state with a different graph (even if they differ by just one 
node). The major advantage of this method is that only vac- 
uum (i.e., offline) squeezing and linear optics are needed, thus 
making this the method of choice in experiments to date lfT9T - 

|23). 

The single-OPO method (U |9) combines all squeezing 
and interferometry into a single optical parametric oscilla- 
tor (OPO), encoding the entire cluster state within a single 
beam. Each mode is an individual frequency within an optical 
frequency comb. The advantage of this method is its scalabil- 
ity. While the initial implementation is more complex than the 
linear-optics method, once the technology is established 11241 - 
1271 . it is in principle much easier to scale up to cluster states 
that are larger by several orders of magnitude [9|. 

The single-QND-gate method [ 10 1 reintroduces the exper- 
imentally challenging Cz gate. But in this case, only one 
such gate is needed because the modes are encoded tempo- 
rally, each traversing the same optical path (but at differ- 
ent times) and each passing multiple times through the same 
optical hardware implementing the Cz gate. This method 
has the additional advantage that the cluster state is extended 
as needed — simultaneously with measurements implementing 
an algorithm on it — in a manner analogous to repeatedly lay- 
ing down additional track in front of a moving train car (a 
"Wallace and Gromit" approach; see footnote in Ref. [10]). 
Such a method eliminates the need for long-time coherence 
of a large cluster state because only a small piece of the state 
exists at any given time. 

The current proposal is based on key elements from the last 
three. The linear-optics method ifTTl is used to generate states 
produced by the single-OPO method [ 8, 9 ] using the temporal- 
mode encoding of the single-QND-gate method [ 10 1. The cur- 
rent method improves over the single-QND-gate method by 
eliminating the need for inline squeezing (in the Cz gate). 



The key observation that makes this simplification possi- 
ble is that the states produced in the single-OPO method 
are Gaussian projected entangled pair states (GPEPS) (28). 
GPEPS states can be described as a collection of entangled 
pairs linked together by projecting adjacent ends down to a 
lower-dimensional joint subspace (schematically, » o » 
becomes a — a — • after projection on the circled nodes). 

The paper is organized as follows. In Sectionlll] I introduce 
the graphical formalism used for the main derivations. Sec- 
tion[llI|describes a temporal-mode GPEPS-based construction 



of a CV quantum wire, and Section IV follows with the same 



for a square-lattice CV cluster state. SectionMconcludes with 
some discussion, including comments on fault tolerance and 
error correction in the CV one-way QC model. 



II. GRAPHS FOR GAUSSIAN PURE STATES 
A. Basic properties 

To describe the states used in this work, I will use the graph- 
ical formalism that covers all Gaussian pure states in a unique 
and unified fashion (29). This formalism is summarized here 
and simplified for the current presentation. The usual graph- 
ical formalism for CV cluster states [6 30, 31] uses graphs 
with real-valued weights and is limited to representing ideal 
CV cluster states, which are necessarily infinitely squeezed 
and thus unphysical. Using the complex-matrix representation 
of Gaussian states ||32ll , the graphical calculus for Gaussian 
pure states [29 1 generalizes these ideal CV cluster-state graphs 
to complex-weighted graphs that can be used to uniquely rep- 
resent any Gaussian pure stateQ including approximate CV 
cluster states. In addition, other types of CV graphs, such as 
the "H-graphs used in the single-OPO method Il9l l33ll . are in- 
corporated directly into the formalism. Local and two-local 
Gaussian unitary operations, as well as local quadrature mea- 
surements, can be visualized as graph transformations. 

To every iV-mode Gaussian pure state \ipz) (ignoring over- 
all displacement and up to an overall phase) we can uniquely 
assign an undirected, complex-wieghted graph on N nodes 
whose adjacency matrix Z is an TV x N complex symmetric 
matrix with positive-definite imaginary part. That is, 



Z = iV + V , 



(1) 



where U = U 1 and V = V 1 are real, and U > 0. All 
Gaussian pure states have a uniquely associated graph of this 
form, and all graphs satisfying these conditions are uniquely 
associated with a Gaussian pure state |29l . This matrix arises 
naturally (up to an overall factor of i) in the position-space 
wavefunction for the state: 



V'z(q) 



(detU) 1 / 4 



T N/4 



exp 



1 T (U- 



*V)q 



(2) 



where q is a column vector of c-numbers. Henceforth, I will 
unambiguously refer to "the Gaussian graph Z." 

All Gaussian unitary operations can be represented by the 
action of a symplectic matrix 



A B 
C D 



(3) 



on a column vector of Heisenberg-picture operators x = 
(p)> where q = (9i)--->?iv) T , and p = ($i,...,pn) t , 
with x = Sx. The new Gaussian graph Z' associated to the 
resulting Gaussian pure state satisfies 



Z i — =— )■ Z' = (C + DZ)(A + BZ) 



(4) 



which can be interpreted as a (generally rather complicated) 
graph transformation rule. However, since all general Gaus- 
sian unitary operations can be decomposed into a sequence 
of local and 2-local Gaussian unitaries chosen from a fiducial 
set, we can build up complicated graph transformations by re- 
peated application of simpler ones. Ref. |29 1 has more details 
on the general case; we will use a simplified version. 

The Gaussian pure state \ipz) satisfies the following nulli- 
fier relation with respect to its graph Z: 



(p-Zq)|^ z >=0, 



(5) 



where the entries in p and q are to be interpreted as 
Schrodinger-picture operators in this context because the 
state \ipx) is indicated explicitlyr] The vacuum (i.e., ground 
state) on all N modes is represented by Z = il, a completely 
disconnected graph with only self-loops of weight i. It is clear 
that Eq. |5]l holds for this state because pj — iqj — —iy/2aj, 
where cij is the annihilation operator for mode j. In the limit 
that U — >• 0, the Gaussian graph Z limits to a real-weighted 
ideal CV cluster-state graph V, and Eq. |5) becomes the usual 
nullifier relation for these states [6|. For our purposes, we 
would also like to represent the H-graph states of the single- 
OPO method [ 33 1 . These states have V = 0. Given an H- 
graph G, the associated Gaussian graph is 



le 



-2oG 



(6) 



where a > is an overall squeezing parameter. 

"H-graphs are incorporated into this graphical formalism 
through the matrix exponential map, as in Eq. (l6j, which is 
normally a highly nontrivial operation. However, it simplifies 
under certain conditions. First of all, note that when G is bi- 
partite (i.e., its nodes can be separated completely into two 
groups such that no edge links nodes of the same group), it 
can be written as 



G 



G J 
Go 



(7) 



2 The formalism only describes the noise properties of Gaussian pure 
states 1291 . Thus, overall displacements are not represented. From now 
on, I will assume that this caveat is understood. 



3 Notice that this is a literal equality, unlike the relation of similar form, 
p — Aq — > 0, for ideal CV cluster states with graph A, which is strictly 
true only in the limit of infinite squeezing (6). 



When G is additionally self-inverse (i.e., G = I), then Gq is 



square and satisfies G Go 
Eq. (|6]l simplifies to 



GqG 



I. In this case, 



Z = % cosh 2a I — i sinh 2a G 



i cosh 2a I 
-i sinh 2a Gq 



-i sinh 2a G 
i cosh 2a I 



(8) 



The essence of the %-graph construction method is that this 
state is equivalent to an approximate CV cluster state after 
phase shifting one of the two sets of nodes (with respect to 
the graph bipartition) J9] 1551 . Performing a Fourier trans- 
form (i.e., phase shift by — |) on the first half of the nodes 
gives |29l 



rji _ ( i scch 2a I tanh 2a G 
I tanh 2a Gq i sech 2a I 



(9) 



which is indeed an approximate CV cluster state since Z' — ► 
G in the infinite-squezing limit (a — > oo)n The fact that 
the ideal CV cluster state approximated by Z has the same 
graph as the 'H-graph G is a peculiar feature of bipartite, self- 
inverse 'H-graphs (9j [34l ; usually the graphs are very differ- 
ent ED. 



B. H-graph states made without an H-graph 

At this point, we should forget about the "H-graph method 
of construction for Z. It doesn't matter how the state repre- 
sented by Eq. ([8) was created — all we care about is that the 
graph for it is given by Z as shown, and we will use graphs of 
this form even when G is only approximately self-inverse (due 
to imperfections at the boundary of a line or lattice). This is an 
important distinction because the simplification of the expo- 
nential map only holds when G is exactly self-inverse. Stray 
connections and imperfections in the H -graph will generally 
"bleed out" via the exponential and contaminate connections 
all over the Gaussian graph Z, which will end up failing to 
produce the state that we want, Eq. (|8), and without an easy 
way to isolate and eliminate the bad links. If instead we start 
directly from Eq. (jSJ, then we'll find — when we phase shift 
the nodes as above — that imperfections in Z are confined to a 
small neighborhood of the imperfections in G. In fact, in this 
case, we have | 



Z' 

{ 



i sech 2a I 



tanh 2a G„ 



\tanh2aGo zsech2a[cosh 2a I — sinh 2aGoG 



(10) 



which reduces to Eq. |9} when G is self-inverse. 

Consider the case in which only a small collection of nodes 
prevents G from being strictly self-inverse. Because the ma- 
trix Z' in Eq. ( fT0] > involves the square of G (through the 
term GoG ), nodes with incorrect edges in G can only af- 
fect other nodes that are at most two steps away (rather than 
everywhere, in the case of the exponential). To see this, de- 
fine G = G + E, where G is bipartite and exactly self- 
inverse, and E is bipartite (with the same partitioning), real, 
symmetric, and sparse. In addition, both G and E have an 
analogous decomposition to that of G in Eq. (fTJ, and there- 
fore Go = Go + Eo, as well. E is a matrix of just a few 
"errors" that cause G to fail to be perfectly self-inverse. We 
have the following relation: 



GoG — (Go 

= I + E G 



Eo)(Go 

T 



Eo) 1 



GqE 



EoE 



(11) 



Now define Po to be a matrix with the same dimensions as Eo 
that is constructed in the following way: first, place a 1 in 
the diagonal entry (Po)jj if row j of E consists of all zeros 
(i.e., if (Eo)jk = Vfc), and place 0's everywhere else; then, 
remove the all-zero rows of Po- Notice that PoEo = 0, which 
gives 



PqGoG P — PqPo =Ij 



(12) 



where the identity matrix on the right is smaller than that in 
Eq. ( [TTj ). Now further define 



I 
Po 



(13) 



which is also not square, where the identity matrix is the same 
size as Eo, and the blocks of are sized appropriatelylj 

Measurements of <jj delete a node j and all of its links from 
a given Gaussian graph [29], which is equivalent to deleting 
the j th row and column from its adjacency matrix. Conju- 
gating Z' from Eq. ( fT0| ) by P deletes from Z' the rows and 
columns that correspond to the nonzero rows of Eo (and, 
equivalently, the nonzero columns of E ): 



P7'P T — i * secn 2al 
Uanh2aP G 



tanh2aGoPo 



i sech 2a I 



i sech 2a I 
tanh2aPGP T , 



(14) 



where Eq. ( fT2| ) has been used. This is an approximate CV 
cluster state with most (but not all) of the graph G intact; the 
troublesome links have been deleted by deleting nodes that 
they attach to. In this case, all of the deleted nodes are taken 



4 Note that this is an approximate CV cluster state despite the fact that it is 
inequivalent to the analogous approximate CV cluster state made by the 
canonical method |4|, whose Gaussian graph would be Z = ie~ 2r l + G. 
In this expression, r is the single-mode squeezing parameter for all the 
nodes that then pass through a collection of Cz gates in accord with G. 
One major advantage of the unified Gaussian graphical formalism is the 
ability to distinguish between these distinct approximants to the same ideal 
CV cluster state. 



5 Alternatively, P can be defined directly in terms of the Moore-Penrose 
pseudoinverse of Eo, which is denoted EJ. Starting with 

(,0 I - E E+J ' 
which is square, we can form P simply by removing the all-zero rows. 



from the second set of the graph bipartition. When these nodes 
are located near the boundary of a regular lattice (or otherwise 
isolated in a graph with local topology), then deleting these 
nodes (or ones just to the "inside" of them) using q measure- 
ments is all that is needed to eliminate the imperfections. 

The point of the preceding analysis is to justify ignoring 
the irregular boundary conditions in the sections that follow. 
In the single-OPO scheme, these boundary irregularities in G 
could not be ignored because they would contaminate the en- 
tire state through the exponential map. In this case, they can 
be dealt with simply by deleting nearby nodes. This is anal- 
ogous to using scissors to clean up the jagged edge of a torn 
piece of paper. 



C. Simplified graphical formalism 

There are only three graph transformation rules that we will 
need for this work: (1) measurements of q, (2) phase-plane 
inversion (i.e., phase shift by if), and (3) interference through 
a 50:50 beamsplitter. All of the graphs we use will be of the 
form of Eq. ([8j, but we will make the following modifications 
to the general formalism |29 1. First, we will not draw the self- 
loops indicated by the diagonal of Z, but we will be mindful 
of them when we derive the simplified graph rules that follow. 
Second, colors indicate the sign of the edge weight but not its 
actual value, with 



— = ^fiC sinh 2a in the graph Z , 



(15) 



where C > 0. Nevertheless, for the graphs used in this work, 
the following relation often (but not always) holds: 



C = max(di,ei2)~ 



(16) 



where di,2 are, respectively, the degrees (number of neigh- 
bors) of the two nodes linked by the edge in questionF] Within 
the approximate CV cluster state Z' that eventually results af- 
ter the phase shifts, these weights become 

— h-> ±C tanh 2a -^-^ ±C in the graph Z' , (17) 

with a — >• oo indicating the infinite-squeezing limit. The fact 
that only the sign of the weights is indicated visually does 
not take away from the rigor of these results; it is merely for 
visual simplicity. When relevant, the actual edge weights will 
be specified. 

Finally, the bipartite nature of Z will be indicated by color- 
ing the nodes black and white, according to the division into 
two sets. This coloring is only a visual distinction and has 



6 This seemingly strange formula arises because of the need for G in Eq. J8l 
to be self-inverse almost everywhere (see Section |lIB) . The exception to 
Eq. (16) is when q measurements are used to delete nodes and their edges 
from the graph; the remaining edges retain their weights before the dele- 
tion. 
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FIG. 1. Simplified graphical formalism for the Gaussian pure states 
represented in this paper, (a) Full graph for a particular Gaussian 
pure state 1 29 1. (b) Simplified graph for the same state, as used in this 
paper. Self loops are not drawn; colored edges indicate the sign of 
the weights as in Eq. l |15| >, with C = -7= (see text) for all edges since 
each edge connects nodes having exactly 2 neighbors; and colored 
nodes indicate the bipartite splitting, with white nodes receiving a 
Fourier transform to convert this to an approximate CV cluster-state 
with the same graph, but now with weights given by Eq. {T71. 



no physical significance on its ownjj The white nodes will 
eventually be Fourier transformed (i.e., phase shifted by — ~) 
to implement the map Z H> Z', as in Eq. (|9j, but since all 
of the graphs we will discuss satisfy Z ~ Z' ~ G up to 
an affine transformation (and neglecting boundary irregulari- 
ties), we will not need a graph transformation for the Fourier 
transform. In addition, since CV cluster states are used for 
one-way QC by making local measurements only, the Fourier 
transform can be absorbed into the eventual measurement pro- 
cedure and becomes simply a change of basisF] An example 
of this simplification of the Gaussian graph formalism is pre- 
sented in Figure [T] 

Below are the graph transformations used in this work. 
These rules are rigorous for the graphs considered in this 
work and are derived from the general transformation rules for 
Gaussian graphs [29], but they have been adapted and simpli- 
fied for the representation used here. These rules function in 
exactly the same way on both black and white nodes. 

Measurement of q — As discussed above, this operation 
simply deletes the measured node and its links from the Gaus- 
sian graph Z. The measured node is indicated by a small ar- 
row: 




/ 



1 In particular, while different colors were used in Ref. [29] to distinguish 
the meaning of nodes in ideal CV cluster-state graphs from their meaning 
in general Gaussian graphs, both colors of nodes in the present paper corre- 
spond to black nodes in Ref. |29|. All of the graphs used here are simplified 
versions of general Gaussian graphs. 

8 Therefore, when we talk about, for example, measuring a white node of the 
resulting CV cluster state in q, we really mean measuring the same node of 
the original state in F^qF = —p. 



Recall that the colors only indicate the sign of the edge 
weights. The actual values of these weights do not change 
with the deletion, but it is worth emphasizing that they may 
no longer satisfy Eq. ([16j. 

Phase-plane inversion (phase shift by w) — This operation 
is also very simple. All edges attached to the inverted node 
(indicated by highlighting) change color: 




This corresponds to the edge weight changing sign but not 
magnitude. 

50:50 beamsplitter — The 50:50 beamsplitter imple- 
ments S B s ( f ) from Ref. 
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(18) 
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Under the conditions of no link between the nodes being in- 
terfered and equal self-loop weights (which is the case with 
all graphs in this paper but is not displayed in the simplified 
formalism), we obtain the following very simple rules: 



>> 




( 



^ 



Notice that the beamsplitter duplicates a single link, but mul- 
tiplies the magnitude of each by a factor of -4= in the process. 
The arrow points from node 1 to node 2 of the interaction, dis- 
tinguishing the effect of the negative signs in Eq. ( fT8| ). When 
an edge is copied in the direction of the arrow, the same color 
is applied to the new edge. When it is copied in the opposite 
direction, the new edge has the opposite color. Even with their 
simplified presentation, all of these graph transformation rules 
are rigorous for the graphs used in this paper. 



III. GPEPS QUANTUM WIRE 

The single-OPO "quantum wire" — a cluster state with a 
one-dimensional topology — actually has a width of two nodes 
and looks like this J5) : 



-iopo 



X~X~S ' C 



(19) 

Notice that only a portion of the full graph is displayedjj We 
can perform a tt phase shift on every other node on the top row 



9 If the full graph has periodic boundary conditions (or if it is formally infi- 
nite), then G for it is self-inverse [9], and Eq. J9j is satisfied. If, however, 



to put this state into a translationally invariant form (ignoring 
node color, which has no physical significance): 



®v^^ * ^S*^* j § v 



;c = i] 

[C = i] . (20) 



The connection to a GPEPS construction is immediate when 
we notice that this is exactly the result of applying the 50:50 
beamsplitter transformation rules to a collection of two-mode 
squeezed states arranged as follows: 




Notice that the order of application of these beamsplitters does 
not matter; the rules commute and result in the same graph. 
The Gaussian graph Z is the desired output of a temporal- 
mode linear-optics construction method. Once in this form, 
measurements of q can be be made on all upper nodes (fol- 
lowing the implicit phase shift by — ^ on all white nodes; see 



Section II C i, projecting the bottom row of nodes into an or- 
dinary (one -dimensional) approximate CV quantum wire [6|: 



-*qw 



[C=\]. (23) 



This GPEPS-based construction method has a simple and 
natural implementation using single-mode squeezers and lin- 
ear optics. This method uses a temporal-mode encoding 1 10|, 
which gives it all the advantages of the equivalent method us- 
ing a single QND gate (in terms of coherence requirements 
and extensibility) but improves on it because all of the re- 
quired squeezing is offline. The method is illustrated in Fig- 
ure |2] 

A few important notes are in order. First, periodic bound- 
ary conditions are used in the single-OPO construction 1 9 1 to 
strictly enforce the requirement that G 2 = I, which other- 
wise fails at the ends. As shown in Section [TlBl however, in 
the construction presented here (which does not rely on ex- 
ponentiating an "H-graph) these "contaminated" nodes can be 
disconnected from the rest of the cluster state by appropriate 
measurements of q on neighboring nodes, just like they are in 
the single-QND-gate method IfTUl . 

Second, it may not be necessary to do the projection of (the 
upper) half of the nodes in q. This was used simply to prove 
that the cluster state with graph G can serve as a CV quan- 
tum wire. Instead, one might be able to use the entirety of G 
as the quantum wire, with quantum information encoded in an 



the graph has at least one boundary, then only Eq. ifToJ is satisfied. As we 
have already seen in Section [TTB] these boundary effects can be dealt with 
by suitably clipping the graph using q measurements. 




FIG. 2. Temporal-mode GPEPS construction of a CV quantum wire 
using passive squeezing and linear optics. Two single-mode squeez- 
ers Si and S2 generate vacuum p- and g-squeezed pulses of light 
(respectively, as shown) at regular intervals At. These pass through 
a simple 50:50 beamsplitter B\, resulting in a two-mode squeezed 
state. (Red arrows point from the first node to the second in Eq. < | 1 8| > 
for each beamsplitter.) The delay loop in the bottom line delays the 
bottom mode by At, allowing it to match up with the top mode of 
the subsequent pair emerging from Bi, resulting schematically in 
the graph shown in Eq. {2TJ. The second 50:50 beamsplitter B2 im- 
plements sequentially each of the transformations indicated by the 
red arrows, resulting in the final graph of Eq. \22\ . These pulses 
head toward detectors Di and D2, which implement the necessary q- 
measurements (phase shifted as appropriate for the white nodes; see 
SectionlnCj, which are indicated by the arrows in Eq. (23}, as well 
as the adaptive measurement-based quantum algorithm to be imple- 
mented on the one-dimensional CV quantum wire. The adaptiveness 
means that subsequent measurement bases generally must be chosen 
based on previous measurement outcomes. Most measurements will 
involve homodyne detection in a basis that must be calculated and 
updated before the arrival of the next pulse, but the ability to divert 
the beam to an efficient photon counter is also required for universal 
single-mode QC l6l . 



appropriate joint subspace of each two-node pair. It is an open 
problem to determine the optimal measurements to perform in 
this case, but initial results indicate the possibility. There are 
several reasons to think this would be useful. First, the final 
cluster state in the construction above has edges with C = |, 
while the two-mode squeezed states all had C = 1. While 
this has no effect on proofs of universality (since a non-unit 
edge weight simply corresponds to squeezing of the quantum 
information as it passes through the cluster), it could have im- 
portant effects on the efficiency of the scheme with respect to 
the amount of initial entanglement in the two-mode squeezed 
states, and weights close to ±1 in the final CV cluster-state 
graph are preferred [ 33 1 . Also, with a joint subspace between 
vertical node pairs now carrying the quantum information, the 
other (unused) joint subspace might play an error correcting 
role. Error correction and fault tolerance for CV one-way QC 
remains an important open problem [28 1 and will be discussed 
further in Section [V] 



IV. GPEPS SQUARE-LATTICE CLUSTER STATE 

The GPEPS techniques developed above for the CV quan- 
tum wire can be adapted to a two-dimensional square-lattice 
CV cluster state, as well. This additional dimension makes 
the state universal for CV one-way QC [4]. We again start 
with the states created by the single-OPO method 1 8 9 ], which 
have the local topology of a square lattice, but with four phys- 



ical nodes per site of the lattice. Analogous to the quantum 
wire, we can make this graph invariant under translations in ei- 
ther lattice dimension by phase shifting the highlighted nodes 
by 7r (which, again, is a local redefinition of basis only): 




(24) 



After doing so, connections between lattice "macronodes" 
(i.e., groups of four nodes) look like this: 




(25) 



This graph is also bipartite (with white and black nodes, as 
shown) and satisfies G = G G = I everywhere except at a 
boundary. 

The GPEPS construction for this state proceeds in two 
steps. We start with the GPEPS construction for a grid of 
unattached, crisscrossing quantum wires: 



Z = 




[C = 1] . (26) 



We then apply two additional beamsplitter interactions to Z : 




C=\], (27) 



which results in the desired graph, shown in Eq. ( |25) . This 
state can be projected down to an ordinary square lattice by 
measuring q on the top three nodes of each macronode (with 
an implicit phase shift by — f on all white nodes; see Sec- 
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FIG. 3. Two equivalent representations of the GPEPS topology of the 
temporal-mode square-lattice CV cluster state illustrated in Eq. {25} , 
adapted from Figure 3 in Ref. 1 10]. Two-mode squeezed states are ar- 
ranged as in Eq. l |26[ >, while red circles represent the four beamsplitter 
transformations indicated in Eqs. (26) and (27) for each macronode, 
as well as the projection using q measurements on three of the mi- 
cronodes within each, eventually resulting in the graph of Eq. ( [28) , 
The two-dimensional square-lattice graph (top) is formally infinite 
in one dimension but finite in the other. This graph can be re- 
drawn as a multiply threaded infinite line graph (bottom). There are 
M additional threadings of the line graph that pass through macron- 
odes M units apart, resulting in a square lattice with vertical dimen- 
sion M. Note that M must be odd to ensure the graph is bipartite. 
(As shown, M = 5.) The dotted links represent additional edges that 
would make the linear version translationally symmetric and equiv- 
alent to a square lattice on a cylinder with one unit of shear in the 
longitudinal direction. Such a family of graphs would still be univer- 
sal for one-way QC because we can measure q on every Mth node 
to delete it (and its links) and "unfold" the graph into an ordinary 
square lattice with a vertical dimension of (M — 1) |6|. 




FIG. 4. Temporal-mode GPEPS construction of a square-lattice CV 
cluster state using passive squeezing and linear optics. Two copies of 
the quantum-wire setup from Figure[2]are used to generate the lattice. 
The upper one has the ordinary delay of At and corresponds to the 
vertical links in Figure[3](top). The longer delay of M At in the lower 
one gives the second threading of the wire and corresponds to the 
horizontal links in Figure [3] (top). Beamsplitters B$ and B4 imple- 
ment the transformations indicated by red arrows in E q. l|26[ >. (Red 
arrows point from the first node to the second in Eq. {T8bfor each 
beamsplitter.) Following this, the 50:50 beamsplitters B5 and Bq 
implement the transformations indicated by red arrows in Eq. (27) , 
eventually resulting in the state with graph Z from Eq. (25) . The 
outputs head to four detectors, which implement the g-measurements 
(phase shifted as appropriate for the white nodes; see Sectionlncl to 
project the state down to an ordinary square lattice, Eq. (28) , as well 
as the adaptive measurement-based quantum algorithm to be imple- 
mented. The adaptiveness means that subsequent measurement bases 
generally must be chosen based on previous measurement outcomes. 
Most measurements will involve homodyne detection in a basis that 
must be calculated and updated before the arrival of the next pulse, 
but the ability to divert the beam to an efficient photon counter is also 
required for universal QC (6). 



tionfflC 




(28) 



Just like for the quantum wire, this projection is useful for 
proving universality of the resulting state, but it may not 
strictly be necessary. Instead, it may be possible to manip- 
ulate quantum information encoded within a macronode as a 
whole (but still using only local measurements) and/or to use 
the additional connections for error correction. Work in this 
direction is ongoing. 

We can use the techniques from the single-QND-gate 
method [ 10] to design a temporal-mode encoding of this state. 
The way this was done in that scheme was to "roll up" the 



square lattice onto a cylinder with a one-unit shift in the longi- 
tudinal direction. This allows the cluster state to be generated 
by double threading a quantum wire. The original threading 
provides the vertical links, while additional connections at a 
spacing of M nodes apart create the second dimension of the 
lattice — see Figure [3] This corresponds to a simple linear op- 
tics circuit, illustrated in Figure [4] As shown in Section II B 



and as in the case of the quantum wire in Section III the "con- 
taminated" nodes at the very beginning of the lattice can be 
disconnected from the rest of the cluster state by appropriate 
measurements of q on neighboring nodes, just like they are in 
the single-QND-gate method ifTOl . 



V. DISCUSSION 

The methods present here for generating a temporal-mode 
CV cluster state using offline squeezing and linear optics com- 
bine the best aspects of several proposals. The temporal-mode 
encoding of the single-QND-gate method [ 10 1 allows finite 
hardware to be used repeatedly by pulsing the single-mode 
squeezer. This method achieves the same result, allowing op- 
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tical elements to be reused by encoding field modes in pulses 
of finite duration. The advantage of this proposal is that the 
need for inline squeezing (Cz gate) is eliminated in favor of 
offline squeezing only, which is sufficient to produce two- 
mode squeezed states when combined with linear optics. This 
was the advantage of the linear optics method [11]. Finally, an 
extensible design combining the two was arrived at by con- 
sidering the states generated in the single-OPO method [8|. 
These states have special properties — specifically, that their 
graphs are bipartite and self-inverse — that allow 50:50 beam- 
splitters to link the two-mode squeezed pairs together. Delay 
loops and phase shifts are the only other ingredients needed to 
generate the states. 

Computing with the states is possible by measuring in q 
the extra "layers" of connections — which exist solely due to 
the use of macronodes in place of individual nodes, and three 
of the detectors in Figure [4] can be assigned simply to this 
purpose. On the other hand, it may be possible to use all four 
detectors in concert to perform C V one-way QC on modes that 
are distributed over the macronode. The optimal measurement 
scheme for using these states is an open area of research. 

The motivation for this is the need for error correction. It 
has been known since the beginning [4 1 that finite-energy CV 
cluster states are, in some sense, inherently faulty due to their 
finitely squeezed nature. Recent results [28 1 show this fault- 
iness to be persistent in the sense that simple encodings of 
qubits or other low-dimensional systems in a finite number of 
CV modes cannot eliminate the errors caused by finite squeez- 
ing. As emphasized in that paper, this does not mean that CV 
cluster states (made by this method or any other) are unusable 
for CV one-way QC. It just means that there is no "magic pill" 
that will eliminate the faultiness introduced by finite squeez- 
ing. Instead, error correction and fault tolerance must be 
addressed from the very beginning if any CV one-way QC 
scheme is to be scalable. While a fault tolerant threshold [ 1 1 
for CV one-way QC is not yet known to exist, promising re- 
sults for qubit-based cluster states [35] and the known ability 
to use CVs for fault-tolerant QC in other contexts [36 1 inspire 
confidence that such a threshold will eventually be found. 

All possible schemes for QC using CV cluster states will 
run into this problem because the faultiness of the states is 
due to finite energy constraints. The current scheme is no 
exception. Nevertheless, there are several reasons to believe 
this is a significant improvement over other optical schemes 
to date. The main reason is that the temporal-mode encod- 
ing limits the need for long-time coherence of the state since 
more cluster is prepared as previous pieces are consumed by 
the detectors. The main theoretical limitation for the current 
scheme is the long delay loop in Figure HI the length of which 
sets the width M of the lattice. While computation can theo- 
retically increase indefinitely in the horizontal dimension, the 



number of horizontal quantum wires linked together (and thus 
the number of encoded qubits or other systems) will be lim- 
ited by how long this delay loop can be made without los- 
ing the ability to phaselock and modematch the interactions 
at the subsequent beamsplitters. Loss and the finite coher- 
ence length of the lasers also come into play, since at least 
2M pulses from each of the four squeezers must continue to 
exist coherently at all times. Despite these restrictions, the 
simplicity of the setup presented in Figure [4] makes it highly 
appealing for new experimental work. Experiments to date, 
which employ the original linear optics method [11|, are cur- 
rently limited to about four modes |fl9l - E2l . 

The main avenues for new research in this area include, 
first and foremost, the experimental implementation of this 
scheme. Beyond that, there remain important questions about 
how best to adapt to the use of macronodes in place of indi- 
vidual nodes in the cluster. This was alluded to several times 
in the text, and the issue relates both to efficiency of imple- 
mentation and also to error correction, fault tolerance, and 
more efficient use of the squeezing resources available. Possi- 
ble extensions include higher-dimensional graphs, such as the 
3D graphs, which are useful in topological one-way QC using 
qubit-based cluster states and result in very high fault-tolerant 
thresholds [35]. The practical limitations of using only finite 
hardware mean that there is a maximum number of encoded 
qubits (i.e., a maximum M) beyond which the scheme will 
cease to be feasible, if for no other reason than the coherence 
length of the laser is exceeded. At this point, either the quan- 
tum wire setup from Figure|2]or the one for the square lattice 
from Figure [4] will need to be concatenated with another such 
system in order to scale up further. The same also applies 
to very long computations if we want them to be done fault 
tolerantly since longer computations will necessarily require 
larger encodings |28|, which reduces to the problem above. 
The precise limitations and means of concatenation remain 
open areas of research. Despite this, given the current state of 
the art in linear optics, we have good reason to believe that the 
method proposed here will provide an important new avenue 
for experimental implementation of CV one-way QC. 
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